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Damaging thermal stimuli trigger long-lasting variation potentials (VPs) in higher plants. Owing to 
limitations in conventional plant electrophysiological recording techniques, recorded signals are 
composed of signals originating from all of the cells that are connected to an electrode. This limitation 
does not enable detailed spatio-temporal distributions of transmission and electrical activities in plants 
to be visualised. Multi-electrode array (MEA) enables the recording and imaging of dynamic 
spatio-temporal electrical activities in higher plants. Here, we used an 8 X 8 MEA with a polar distance 
of 450 ]im to measure electrical activities from numerous cells simultaneously. The mapping of the data 
that were recorded from the MEA revealed the transfer mode of the thermally induced VPs in the leaves 
of Helianthus annuus L. seedlings in situ. These results suggest that MEA can enable recordings with 
high spatio-temporal resolution that facilitate the determination of the bioelectrical response mode of 
higher plants under stress. 



Electrical signals have been found in many plant species, including lower plants and both sensitive and 
non-sensitive higher plants. These signals are involved in many plant biological processes such as 
respiration', water uptake^, leaf movement"*, and biotic stress responses*. The two main electrical signals 
that are found in plants are action potentials (APs) and variation potentials (VPs). While APs are rapid 
signals that are induced by non-damaging stimuli\ VPs are slower signals that are induced by damaging 
stimuli*"' (e.g., burning and cutting). Indeed, VPs are sometimes called slow-wave potentials because their 
transmission speed is typically slower than that of APs". Additionally, VPs decrease in magnitude as they 
spread away from stimulated sites'. Both types of signals carry information to other cells, tissues, and organs, 
where they trigger appropriate responses. However, the modulation of VP amplitude, which is accompanied 
by interdependent ion fluxes'", increased ethylene emissions", and induced proteinase-inhibitor gene 
expression'"', may enable the VPs to carry more information than the APs. In fact, VPs that are induced 
by different stimuli have distinguishing signal characteristics and propagation patterns'^. However, compared 
with APs, the patterns and pathways of VP transmission are less well known. Because of limitations in 
recording wound-induced VPs by surface electrodes, distinguishing changes in electrical potentials from 
multiple cells or tissues is impossible in plants, as is the mapping of VP propagation. 

Multi- electrode arrays (MEAs) can overcome these limitations and enable the monitoring of bioelectrical 
activities at multiple sites with high spatio-temporal resolution'''. MEAs have been used in vitro for more than 40 
years to measure the electrical activities of heart cells from cultured embryonic chicks"". A MEA device consists of 
an array of microelectrodes (typically 10-50 |im in diameter) that are connected to an external circuitry that 
enables the recording of electrical activities from excitable cells (usually neurons or cardiomyocytes) and tissues. 
In recent years, significant progress has been made in this technique, and highly reliable MEAs have been used in 
the extracellular and intracellular recording of animal cells """ and maize root apices™. Here, we employed an 8 X 
8 MEA to measure VPs in H. annuus leaves in vivo and in situ and investigated the propagation of VPs in a region 
of interest (ROI). 



SCIENTIFIC REPORTS | 4:5435 | DOI: 1 0. 1 038/srep05435 



1 




) ; , , , , , , ; , , 1 1 

0 10 20 30 40 50 60 70 80 90 100 110 120 
Time (s) 



Figure 1 | Coupling of a H. annuus leaf and electrodes of the MED64 array, (a) The leaf was placed on the electrode panel, and a U-shaped platinum plate 
with a grid of nylon threads was placed on the top of the leaf A 25-g weight that was sealed with a plastic membrane was then placed on the platinum plate. 
A 10 mM KCl solution was placed in the electrode chamber, (b) A microscopic view of the planar electrodes and leaf after recording, showing the 
close connections between them; the black dots are the electrodes, (c) Baseline activities of 64 electrodes after coupling, (d) Overlapped baseline signals of 
the 64 channels that are shown in (c), indicating that the noise level was less than ±17 |iV. 



Results 

Noise level of the MED64 array. Figure la shows the coupling of 
electrodes in the array and a H. annuus leaf. Figure lb shows a close- 
up view of the connection between the planar electrodes and the leaf. 
In the present study, we found that the ion concentration of a 
solution in the probe chamber had an effect on the noise level of 
the recording results. The recording results with different KCl con- 
centrations (1 mM, 2 mM, 5 mM, and 10 mM) and the standard 
buffer solution (2 mM KCl, 1 mM CaClz, 1 mM MgClz, 50 mM 
mannitol, and 2.5 mM MES/NaOH buffer, pH 5.7) demonstrated 
a noise level of approximately ± 50 |J,V in 1 mM KCl, approximately 
±17 |iV in 10 mM KCl, and approximately ±18 |iV using a stand- 
ard buffer solution (see Supplementary Fig. SI). Thus, 10 mM KCl 
and the standard buffer solution were used in the experiments, 
respectively. Figure Ic and Figure Id show the baseline and noise 
level of the MED64 array. The noise level was less than ±17 |J,V at a 
band-pass range of 0.1-1 KHz and a sampling frequency of 20 KHz 
at 10 mM K+. 

Recording of variation potentials using MEA and analyses of 
signal transmissions. We recorded heat-induced electrical respon- 
ses from 45 different plants in distinct experiments (38 plants for 
recording in 10 mM KCl and 7 plants for recording in standard 
buffer solution) and observed similar transmission patterns for all 
of the experiments. A typical example of a response to thermal 
stimulation is illustrated in Fig. 2. Figure 2a indicates the relative 
position between the planar microelectrodes and the H. annuus 
leaf We stimulated the right leaf and recorded the left leaf 
Figure 2b shows the field potential that was induced by thermal 
stimulation, and nearly all 64 of the electrodes successfully 



recorded VPs. In Fig. 2c, the signals from Fig. 2b are overlapped, 
showing that the VPs occurred with a mean latency of 55 ± 5 s after 
stimulation onset, which lasted 30 ± 2 s. Similar response times have 
been reported in other studies'". The signal amplitudes were in the 
range of 50-400 |.iV Vp-p. Clear delays can be seen among the 
different VP traces, indicating that VP transmission patterns 
differed depending on the ROI. To determine the propagation 
directions and transmission routes of the VPs within the recording 
region, we spatially mapped each VP according to its placement 
within the MEA of the electrode that recorded it. The direction of 
propagation is shown by the colour coding of the VPs according to 
their onset time in Fig. 2d (blue, early onset; yellow/green, mid-range 
onset; red, late onset). The colour coding revealed that the VPs first 
occurred on the right side (electrodes 8, 16, 24, 32, 40, 48, and 56) of 
the MEA and gradually progressed to the upper left side (electrode 1 ) . 
The time delay between the first and last VP onset (electrodes 8 and 1, 
respectively) was 6 s. The propagation velocities of the VPs in any 
direction can be easily deduced from the time delay between the 
detection of the signal at the two selected electrodes (shown by 
arrows in Fig. 2e). The propagation speeds ranged from 0.7 mm/s 
to 2.5 mm/s and occurred in different directions. The VP that was 
transmitted along the main propagation route (red arrow. Fig. 2e) 
had a slower propagation velocity than those of the VPs that were 
transmitted along the branching routes (blue arrows. Fig. 2e). 
Considering the relative positions of the electrodes on the leaf 
(Fig. 2a), the VPs were clearly induced in the right leaf, trans- 
mitted along the petiole to the left leaf, and subsequently travelled 
across the surface of the left leaf Figure 2f uses time-lapse pseudo- 
colour images to illustrate the amplitude changes that occurred as the 
VPs propagated across the leaf The potential change initially 
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Figure 2 | Spatio-temporal propagation of waves in a leaf neighbouring the heat-stimulated leaf, (a) Relative position of planar microelectrodes 
(left leaf) and location of stimulation (right leaf), (b) After thermal stimulation, nearly all 64 of the electrodes recorded VPs. The white arrow indicates an 
artifact that was observed during the thermal stimulation, (c) Overlapping signals from (b). The VP onset was delayed approximately 50 s after 
stimulation, and the VP duration was approximately 30 s. (d) The VP propagation was visualised through the colour-coding of the VP onset at each 
electrode. The VPs reached the electrodes that are depicted in blue first and then propagated to the yellow regions and finally to the red regions, 
(e) Average VP propagation velocities in the different directions. The velocity equals the distance between two numbered planar electrodes on the arrow 
divided by the time delay between signal detection at the two electrodes, (f) Maps of the voltage distributions for the propagating waves at individual 
snapshots in time. Changes in voltage partially match the propagation direction (right to left) that is shown in (d). 



occurred in the right recording side (blue squares, t = 64.5 s) and 
then spread to the left, which is a pattern that partially matches the 
transmission direction that is seen in Fig. 2d. 

The results from the experiment shown in Fig. 3 cross-validate 
those from the experiment shown in Fig. 2. Because thermal stimu- 
lation damages the plant, we used different H. annuus samples in this 
experiment. Compared with Fig. 2a, the relative positions for record- 
ing and stimulation were opposite. Again, nearly all 64 of the elec- 
trodes recorded VPs that were induced in the left leaf by the thermal 
stimulation of the right leaf (Fig. 3b), and the delay time was approxi- 
mately 40 s (Fig. 3c), which is similar to the results obtained in the 
first experiment. Although the 20-s duration was shorter than that in 
the first experiment, the times were comparable. Additionally, varia- 
tions in onset delay can be seen across the electrodes. For example, 
the delay time was 5 s between electrodes 25 and 40 (Fig. 3d). The 
propagation speeds were also similar to those that are illustrated in 
Fig. 2, measuring from 0.8 mm/s along the main route (red arrow in 



Fig. 3e) to 2 mm/s along a branching route (blue arrow in Fig. 3e). 
Figure 3d and Figure 3f show that the VPs first reached electrodes on 
the left side of the grid and then travelled to the right side, which was 
again consistent with the experimental setup, indicating that the VPs 
spread with characteristic patterns across the leaf surfaces and from 
one leaf to another. Additionally, we recorded the leaf VPs in exactly 
the same way for the MED probe with the standard buffer solution. 
The results are consistent with those that are shown in Fig. 2 and 
Fig. 3 (see Supplementary Fig. S2). 

Furthermore, we used conventional inhibitors of H*-ATPase and 
ion channels, including sodium orthovanadate. La'*, and tetraethy- 
lammonium ion (TEA*), to examine this MEA recording method. 
The VPs of the leaves that were treated with inhibitors as recorded by 
MED64 are shown in Supplementary Fig. S3. The profiles in 
Supplementary Figs. S3a, b, and c show the VP recorded approxi- 
mately 70 s after heat stimulation. The maximum delay times of the 
leaves that were treated with sodium orthovanadate increased to 7.3 
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Figure 3 | Spatio-temporal propagation of waves in a leaf neighbouring the heat-stimulated leaf; cross validation of results using a different plant. 

(a) Relative positions of the planar microelectrodes (right leaf) and the stimulation location (left leaf), (b) After thermal stimulation, nearly all 64 of the 
electrodes recorded VPs. The white arrow indicates the artifact that was observed during the thermal stimulation, (c) Overlapped signals from (b). 
(d) The VP propagation was visualised through the colour coding of the VP onset at each electrode. The VPs reached the electrodes that are depicted in 
blue first, then propagated to the yellow regions, and finally moved to the red regions, (e) Average VP propagation velocities in different directions. The 
velocities were calculated as described in the caption of Fig. 2e. The propagation velocity was slower along the main route than along the branching routes, 
(f) Maps of the voltage distributions for the propagating waves at individual time points. Changes in voltage match the propagation direction 
(left to right) that is shown in (d). 



± 1.5 s, 6.3 ± 1.4 s (La'+), and 7.2 ± 1.3 s (TEA+), compared with 
the control delay time of 5.5 ± 1.7 s. Although there is a slight 
difference between the maximum delay time in Supplementary Fig. 
S3d and that in Fig. 2, the mean velocities ranged from 0.6 mm/s in 
the main route to 2 mm/s in the branch route, which are smaller than 
those shown in Fig. 2 and Fig. 3. The maps of the voltage distributions 
for the propagating waves at individual time points are displayed in 
Supplementary Fig. S3f, which shows that the amplitudes of the VPs 
after inhibition are much smaller than those of the controls. All of the 
blockers were shown to significantly reduce the amplitudes of the 
VPs compared with the control (p < 0.05) (Fig. 4a) and increase the 
maximum delay times among the signals, although these increases 
were not significant (Fig. 4b). As shown in Fig. 4, the H^-ATPase and 
ion channel blockers altered the amplitudes of the VPs from control 
values of 249 ± 51 |j,Vto71 ± 13 |,iV (1 mM sodium orthovana- 



date),67± 15 |iV(l mM La'+), and 56 ± 15 |xV(10 mMTEA+).In 
addition, the recorded leaf petioles were killed by the 10% para- 
formaldehyde solution, no VPs were recorded for those leaves 
(n = 8) (see Supplementary Fig. S4). 

Discussion 

Using a MEA, we recorded the VPs in the leaves of H. annuus in vivo 
at multiple sites simultaneously. In contrast to surface electrodes, 
MEA has high spatio-temporal resolution that can detect changes 
in potentials in a 3.55 mm X 3.55 mm area using 64 planar 
microelectrodes. 

In our experiments, the VP amplitudes were in the range of 50- 
400 |iV. Compared with traditional extracellular recording using 
surface electrodes, in which VP amplitudes usually range from sev- 
eral millivolts to tens of millivolts'"'^', the amplitudes that were 
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Figure 4 | The effects of pharmacological blockers on the amplitudes and propagations of the VPs that were recorded using MED64. (a) In the control 
group, no blockers were added to the leaf, and the mean amplitude ofVP from the 64 multi-electrodes was 236 ± 68 |^V. After the leaf was treated using 
the plasma membrane H*-ATPase inhibitor sodium orthovanadate, the Ca^* channel inhibitor La"**, or the K* channel inhibitor TEA*, the mean VP 
amplitude was reduced to 71 ± 13 ^V, 67 ± 15 )iV, and 56 ± 15 ^V, respectively. Data are presented as the mean ± S.E.M. unless stated otherwise 
(n = 6 per group). *p < 0.05 (Student's t-test) for the treatment groups compared with the control, (b) There were also increases in the maximum time 
delays of the signals (1 mM sodium orthovanadate, 7.3 ± 1.5 s; 1 mM La'*, 6.3 ± 1.4 s; 10 mM TEA*, 7.2 ± 1.3 s) compared with the control delay time 
(5.5 ± 1.7 s). Data are presented as the mean ± S.E.M. The maximum time delay between the treatment groups and the control was not significant. 



observed here were much lower. The differences in VP amplitudes 
between the MED64 system and traditional methods are a result of 
the tissue/microelectrode interface model that is used in the planar 
microelectrode array but not the traditional methods^^. Specifically, 
the position of the reference electrode affects the magnitudes of the 
VP amplitudes that are recorded^'. The recording electrodes and 
reference electrodes of the MED64 probe are arranged in the same 
planar glass plate^'*, which records the voltage that is generated 
between the leaf contacting the electrodes and the bath solution 
(10 mM KCl) or the standard buffer solution. These results show 
lower amplitudes than those that were recorded using traditional 
methods. Although many factors and parameters may affect the 
performance of the system, the most important factor for the accu- 
racy of MEA is to ensure that the system has an acceptable signal/ 
noise ratio^*'^^. Although the VP amplitudes were only between 
50 \xY and 400 |iV, we were still able to distinguish the details of 
the VPs because the noise level of the MED64 system was extremely 
low (±17 |tV). We found that there was no effect on the VPs of the 
leaves and VP propagations using the MED probe with 10 mM KCl 
(n = 38) or a standard bath solution (n = 7). Thus, the results are 
reliable and comparable following non-invasive and nondestructive 
in situ recording processes using different solutions. Incidentally, 
Furch et al""" noted that a high concentrations of external KCl 
(100 mM) evoked a rapid AP-like response; thus, it is important to 
prevent any breakage from the measured leaf from occurring. 

In general, the major long transmission routes of the VPs can be 
determined using traditional electrophysiology methods, but 
because of limitations in the numbers and sizes of electrodes, there 
are few studies examining VP transmission routes at the cellular 
level, and visualising the transmission of electrical signals using tra- 
ditional methods is not possible. In the present study, we found that 
the use of MEA may solve these problems. Our results (Figs. 2d-e 
and Figs. 3d-e) revealed that the transmission routes of VPs in H. 
annuus leaves did not follow a line or travel from point to point but 
spread across the surface in multiple directions. The data that are 
shown in Fig. 2e and Fig. Be suggest the presence of two types of 
transmission routes (the blue and red lines) with obvious differences 
in transmission velocities. Signals travel slowly along the main vein 
(parallel to the leaf) and much faster along the smaller veins that 
branch off towards the edges of the leaf However, signals do not 



seem to travel directly between the side branches; they arrive at each 
part of the leaf via the main vein and then follow the branches 
outward. Considering that different tissues (e.g., veins, epidermis, 
guard cells, and mesophyll) were covered by the electrode array, it 
is likely that the VPs spread at different velocities depending on the 
leaf tissue. Our findings also support the idea that VPs are primarily 
generated by wounding, which induces a hydraulic wave that is 
transmitted through the xylem, causing a local electrical response 
in neighbouring symplastic cells that propagates faster than the res- 
ponse that occurs in the main vein. The signal can be transmitted 
over short distances from cell to cell through the plasmodesmata, and 
after it reaches the phloem, it can also propagate over long distances 
along the sieve-tube plasma membrane"*'^''^". 

The possibilities and advantages of visualising VP transmission in 
plants can also be seen in Fig. 2d and Fig. 3d. VP propagation from 
one side of the detectable region to the other was indicated by the 
analysed data from the MEA. Complex neuronal networks found in 
animals do not exist in plants, and many signal patterns, such as 
bursting or synchronisation that have been widely studied in animal 
nerves''^"'' and brains'", have not been discovered or investigated in 
plants in situ. Considering the variations in the features of electrical 
signals in plants, different information (perhaps corresponding to 
different environmental stimuli, such as cold, heat, or cutting) might 
be transmitted through different VP signal types. Wounding- 
induced VPs have been associated with many responses in plants, 
including gene expression in tomato"''^ and Arabidopsis thaliana'^, 
the cessation of mass flow and callose deposition in Cucurbita max- 
ima?*. We also observed callose deposition in the present study (see 
Supplementary Fig. S5). In addition to VPs, the electrical responses 
to cooling stimulation were also measured using MEA (see 
Supplementary Fig. S6). Because a MEA can detect electrical signals 
at multiple sites simultaneously, this technique holds great potential 
for identifying different signal transformations. Furthermore, when 
combined with fluorescent dyes or other techniques, it may be 
applied to study the relationship between electrical activity and 
physiological responses at multiple sites in the plant. 

In many higher plants, VPs have been recorded in leaves^', stems''\ 
sieve tubes'", and petioles'' using surface electrodes, aphid stylets, 
and microelectrodes. Based on these investigations, the major 
notions regarding VP transmission are focused on the hydraulic 
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wave^^ and chemical-inducing^**'^^ theories, in which the xylem and 
phloem are suggested to be the long transmission routes for VPs. In 
the present study, we used conventional inhibitors to investigate the 
factors influencing VPs and VP propagation. Our results agree with 
those that have been reported in the literature**". In particular, our 
observations are consistent with reports describing effects of Ca^"^ 
channel blockers^^ on electrical signals. It has been proposed that an 
electrical potential wave (EPW)-associated Ca^"^ influx may trigger 
the generation and amplification of countless long-distance sig- 
nals^^'^^, and these results have encouraged exploration of the ion 
mechanisms that are involved in electrical activity in the phloem. 
Moreover, the mapping of VP propagation with high spatio- 
temporal resolution enables us to determine the distribution features 
of the VPs and VP propagation patterns. Given the strong research 
interest in vPs^^'^^^'^^"**^, this method is a potentially useful tool 
for the elucidation of the electrical signals that occur following 
wounding. 

In this study, VPs in H. annuus leaves were recorded at multiple 
sites simultaneously using a multi- electrode array. Our results 
revealed that VPs have a complex surface-transmission pattern in 
H. annuus leaves and indicated the potential of using MEA to visua- 
lise VP transmission in plants. Our research suggests that MEA can 
provide high spatio-temporal resolution recordings to determine the 
bioelectrical responses of higher plants under stress. 

Methods 

Plant materials. The H. annuus plants that were used in this study were grown in a 
greenhouse under a 16 h light/8 h dark photoperiod regime at temperatures of 30"C 
and 25 'C during the day and night, respectively. The growth medium was comprised 
of 70% vermiculite and 30% humus and contained all of the essential nutrients. The 
plants that were used for the experiments were 2-3 weeks old. 

MEA system. We used an MED64 {Alpha MED Scientific Inc., Ibaraki Osaka, Japan) 
MEA system. This system has 64 planar microelectrodes that are arranged in an 8 X 8 
array in the centre with four reference electrodes that are symmetrically situated at the 
four corners of the probe plate (MED-P545A, Alpha MED Scientific Inc.). Each 
microelectrode size is 50 |j.m X 50 |im, and the distance between adjacent electrodes 
is 450 |j.m. Thus, the array can cover an area of approxunately 3.55 mm X 3.55 mm. 
According to the MED64 Hand Book form Alpha MED Scientific Inc. (Japan), the 
four reference electrode lines flow together and then converge on a common ground 
in the MED64 amplifier. Additionally, the user manual suggests that the platinum 
wire ring functions as an additional reference electrode to suppress noise and stimulus 
artifacts. A glass ring (I.D. — 22 mm, O.D. — 25 mm, height — 10 mm) was attached 
at the centre of the plate to form a recording chamber. An additional platinum ring 
(O.D. — 21 mm) reference electrode was placed at the bottom of the recording 
chamber to connect with the ground, but it did not contact the leaves during 
recording. The sampling frequency and bandwidth were set to 20 kHz and 0.1 Hz to 
100 kHz, respectively. The input impedance of the amplifier was 100 MQ. 

Preparation of multi-electrode array. Procedures for the preparation of the 
introduced multi-electrode probe (MED-P545A, Alpha MED Scientific Inc.) were 
previously described by Zhao et al*^. The MED64 probe was treated with 0.1% 
polyethyleneimine (Sigma) in 25 mM borate buffer (pH 8.4) overnight at room 
temperature*^. Immediately before use, the probe surface was rinsed three times with 
sterile distilled water. 

Coupling of if. annuus leaves and MED64 planar microelectrodes. For recording, a 
leaf was fixed in situ in a recording chamber and placed directly on the multi- 
electrode array, which was connected to the 64-channel system using the MED 
connector. The chamber was filled with 1 mL of a 10 mM KCl solution or standard 
buffer solution, which immersed the plate electrodes and platinum ring. The recorded 
part of the leaf was also immersed in the solution. To achieve good contact between 
the electrodes and the leaf, a U-shaped platinum plate with a grid of nylon threads was 
placed on the top of the leaf, and a 25-g weight that was sealed with a plastic 
membrane was placed on the platinum plate. 

Inhibitor Treatment. Inhibitors were added separately to the experimental chamber 
containing a standard buffer solution of 2 mM KCl, 1 mM CaCl2, 1 mM MgCl2, 
50 mM mannitol, and 2.5 mM MES/NaOH, pH 5.7. The concentrations of inhibitors 
were as follows: Ca^"^ channel inhibitor, 1 mM LaCls; K"^ channel inhibitor, 10 mM 
tetraethylammonium chloride; and H^-ATPase inhibitor, 1 mM sodium 
orthovanadate. Prior to heat stimulation, the duration of incubation for the leaf in the 
standard solution containing the inhibitors was at least 1 h. After the treatment, all of 
the recording conditions that were used were the same as those for the untreated 
control plants. 



Stimulation. Thermal stimulation was applied to a nearby leaf, which was separate 
from the recording leaf Stimulation was performed by passing a flaming match under 
the leaf tip for 3 s in an area measuring approximately 1.5 cm^. The distance between 
the stimulation point and recording region was 8-10 cm. 

Data analysis. Data were collected using MED64 Mobius acquisition software 
(Witwerx Inc., Tustun, CA, USA) and analysed and visualised using the Matlab 
software (Matlab R2009). Data are presented as the mean ± S.E.M. A Student's t-test 
was used to analyse the differences among groups, andp < 0.05 was considered 
significant. 
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